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Abstract 

In this paper we consider FRW cosmology in modified gravity which contain arbitrary functions 
f(<fi). It is shown that the bouncing solution appears in the model whereas the equation of state 
(EoS) parameter crosses the phantom divider. The reconstruction of the model is also investigated 
with the aim to reconstruct the arbitrary functions and variables of the model. 

PACS numbers: 04.50.Kd; 98.80.-k 

Keywords: Modified gravity; Bouncing universe; oj crossing; Stability condition; Reconstructing. 



*Elcctronic address: hosseinf@guila n7ac.ir| 
^Electronic address: fmilani@guilan.ac.ir 



1 



1. INTRODUCTION 



There are many cosmological observations, such as Super-Nova la (SNIa) Wilkinson 
Microwave Anisotropy Probe (WMAP) Q, Sloan Digital Sky Survey (SDSS) fl, Chan- 
dra X-ray Observatory 4j etc., that reveal some cross-checked information of our universe. 
They suggest that the universe is spatially flat, and consists of approximately 70% dark 
energy (DE) with negative pressure, 30% dust matter (cold dark matters plus baryons), and 
negligible radiation, and also the universe is undergoing an accelerated expansion. 

Recent observations have determined basic cosmological parameters in high-precisions, 
but at the same time they posed a serious problem about the origin of DE. The combined 
analysis of SNIa js]], that is based upon the background expansion history of the universe 
around the redshift z < 0(1), galaxy clusters measurements and WMAP data, provides an 
evidence for the accelerated cosmic expansion jg]. The cosmological acceleration strongly 
indicates that the present day universe is dominated by smoothly distributed slowly varying 
DE component. The constraint obtained from SNIa so far has a degeneracy in the EoS 
of DE 7|. To many people's frustration, the ACDM model with an EoS u — — 1 has 
been continuously favored from observations. This degeneracy has been present even by 
adding other constraints coming from Cosmic Microwave Background (CMB) and Baryon 
Acoustic Oscillations (BAO) 91. The modern constraints on the EoS parameter are around 

nn 

the cosmological constant value, uj = — 1 ± 0.1 [6|]-[10|] and a possibility that u is varied in 
time is not excluded. From the theoretical point of view there are three essentially different 
cases: u > — 1 (quintessence), u = —1 (cosmological constant) and u < — 1 (phantom) 
( 111-I14J and refs. therein). 

nn 

The models of DE can be broadly classified into two classes [13, [16| . The first corresponds 
to introducing a specific matter that leads to an accelerated expansion. Most of scalar field 



models such as quintessence 



17| and k-essence 



belong to this class. The second class, 



that in this paper we consider, corresponds to the so-called modified grav ity models such as 
f(R) gravity 19| . scalar-tensor theories 20[] and brane-world models 2l|]. In order to break 
the degeneracy of observational constraints on u and to discriminate between a DE models, 
it is important to find additional information other than the background expansion history 



of the Universe 



22]. 



In second classification, modified gravity 23j suggests fine alternative for DE origin. 
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Indeed, it may be naturally expected that gravitational action contains some extra terms 
which became relevant recently with the significant decrease of the universe curvature. The 
modified gravity can be obtained in two ways, first by replacing scalar curvature R, or with 
f(R), in the action which is well known as modified gravity ,or f(R) modified gravity, and 
second by considering additional curvature invariant terms like Gauss-Bonnet (GB) term. 
Another modification of GR is modified Gauss-Bonnet gravity 24| which is obtained by 
inserting a function of GB invariant f{G), in the Einstein-Hilbert action. A number of 



which explain the 



metric formulation of modified f(R) gravities has been proposed 
origin of cosmic acceleration. Particular attention is paid to f{R) models 29]- 32] with the 
effective cosmological constant phase because such theories may easily reproduce the well- 
known ACDM cosmology. Such models subclass which does not violate Solar System tests 
represents the real alternative for standard General Relativity. 
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On the other hand, the Friedman equation forms the starting point for almost all in- 
vestigations in cosmology. Over the past few years possible corrections to the Friedman 
equation have been derived or proposed in a number of different contexts, generally inspired 
by brane- world investigation 3JL |35fl . These modification are often of a form that involves 



the total energy density p. In 



36] , multi-scalar coupled to gravity is studied in the context 



of conventional Friedman cosmology. It is found that the cosmological trajectories can be 
viewed as geodesic motion in an augmented target space. 

There are several phenomenological models describing the crossing of the cosmological 



constant barrier 
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38] . Most of them use more then one scalar field or use a non- minimal 



coupling with the gravity, or modified gravity, however we use both of them in our paper. 
In two-field models one of these two fields is a phantom, other one is a usual field and the 
interaction is non-polynomial in general. It is important to find a model which follows from 
the fundamental principles and describes a crossing of the u = — 1 barrier. 

A bouncing universe which provides a possible solution to the Big Bang singularity prob- 
lem in standard cosm olog y has recently attracted a lot of interest in the field of string theory 



and modified gravity 39|, |40[. In bouncing cosmology, within the framework of the standard 



FRW cosmology the null energy condition (NEC) for a period of time around the bouncing 
point is violated . Moreover, after the bouncing when the universe enters into the hot Big 



Bang era, the EoS parameter cu in the universe transits from u < — 1 to > — 1 



4l|- 



In this paper, in section 2 we study the dynamics of the FRW cosmology in modified 
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(non-local) gravity. We discuss analytically and numerically a detailed examination of the 
conditions for having u across over —1. The necessary conditions required for a successful 
bounce is discussed in this section as well. Section 3 describes our model reconstruction and 
present the results. Finally, we summaries our paper in section 4. 

2. THE MODEL 

We start with the action of the non-local gravity as a simple modified garvity given by 

r ( m 2 i 

S = j d'x^-g \-fR{i + f{n- x R)) j- (l) 

where M p is Plank mass, / is some function and □ is d'Almbertian for scalar field. Generally 
speaking, such non-local effective action, derived from string theory, may be induced by 
quantum effects. A Bi- scalar reformation of non-local action can be presented by introducing 
two scalar fields <p and ip, where changes the above action to a local from: 

r _ ( M 2 1 

S = j ( i 4 x v ^|^[ J R(l + /(0)) + ^(n0- J R)]j, (2) 

where ip, at this stage, plays role of a lagrange multiplier. One might further rewrite the 
above action as 

/( M 2 1 
d'x^-g }^ [R{1 + m -^- d^4>} j , • (3) 

which now is equivalent to a local model with two extra degrees of freedom. By the variation 
over ip, we obtain Ucf> = R or = D" 1 /?, where / is a function of scalar field in the model. 

Now in a FRW cosmological model with invariance of the action under changing fields 
and vanishing variations at the boundary, the equations of motion for only time dependent 
scalar fields, <fi and ip, become 

4> + 3H<p + R = 0, (4) 
+ 3Hijj - Rf = 0, (5) 

where R = 12H 2 + 6H, H is Hubble parameter and /' = • Variation of action ([3]) with 
respect to the metric tensor g^ v gives, 

= l -g„ v {R(l + f-ip)- d p m} - R, u (l + / - V) + \{d^d u <j> + d^d^) 

- Ovd-v^xj-vO- (6) 
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The 00 and ii components of the equation ([6]) are 



-3H 2 (1 + f - + \^ - 3H{fj> - </>), 



= (2H + 3H 2 )(l + f-^) + ^<P 
Equations © and (jSJ) can be rewritten as 

3H 

2H + 3H 2 



(7) 
(8) 



A#-3if(/'0-V) 



(1 + /-^) 
|^+(^ + 2i/|)(/-^) 



(9) 
(10) 



i + f-4') 

2M p 2 ' Ult! 

right hand side of the equations (Q and (flOl) can be treated as the effective energy density 
and pressure: 



Comparison with the standard Friedman equations H 2 = |^r, and 2H + 3H 2 — ^jp, 



M 2 



\i>j>-3H(f'i>-i>) 



(1 + /-V0 ' 



M 2 

p 



(11) 
(12) 



+ 

Using Eqs. (jlj) and (jSJ) and doing some algebraic calculation we can read the effective energy 
density and pressure from the above as, 



Peff 



M 2 r i 

J-^-3F(/V-V0 



1 + 12 



(13) 



M 2 



Pe// 



1 + / - V - 6/' 2 



7<^0 + /> 2 -^(/V-^) + 



1 + 



Now by using Eqs. (JT3|) and (jHJ) the conservation equation can be obtained as, 

p ef f + 3H p eff (1 +u) = 0, 
where a> = is the EoS parameter of the model. Also from Eq. (J7]) we obtain, 



H _ -(M-ip) / 6 #(l + /-^) 

"2(1 + / --0)1 V 9 W-^) 2 



(14) 



(15) 



(16) 



At this stage we study the cosmological evolution of EoS parameter, u , and show that 
analytically and numerically there are conditions that cause the EoS parameter crosses the 
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phantom divide line (u — > — 1). Let's see under what conditions the system will be able to 
cross the barrier of uj = — 1. In order to do that, one requires Peff + Peff to vanish at a point 
of (0o, "00 ) an d change the sign after the crossing. This can only be achieved by requiring 
H((fi ,if) ) = and H has different signs before and after the crossing. 

To explore this possibility, we have to check the condition J^(p e // + Peff) 7^ when 
cu — > — 1. Using Eqs. (|T3|) and (fl4|) in second Friedman equation, Eq. ( flT)]) gives, 



■ _ 24fH 2 + mfj - erg + 

2(1 + / - ^ - 6/') 1 j 

Also we have ^§(p e // + Peff) = —1M 2 p E or, 

0(4#/' - 2/"0 - V>) - ^(4# + 0) + <j>(4Hf"(6H + 0) - f"(j) 



H 



2(1 + / - V - 6/') 
, * f 24/'/r + 0(/' + 6/")-^A 

+ ^ l + f-t-er )■ (18) 

One can find that in order to have w-crossing one of the following conditions might be 
satisfied when u — > — 1 and H ^ 0. 

• (a) -0 = and 0^0 

• (b) = and ip ^ 

• (c) = and -0 = 

In the first case, (a), we have 

and, 

20(2tff - f 0) - $(4# + 0) + (4Hf"(6H + 0) - f "0 2 ) 
H = 2(l + /-V-6/') (21) 

Therefore the conditions for having u across over —1 are: (a-1) ^ and 2///' 7^ /"0 
when other terms can be neglected, (a-2) ip 7^ however the first and third terms can be 
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vanished, (a-3) 2<j>{2H f - f" 4>) ^ -0(4#+0) where third term is zero or (a-4) AH f"(6H+<j>) ^ 
f"'(j) 2 or if = then /" 7^ and if /" = then 7^ 0, when the first and second terms 
can be vanished in addition to the f"(p 2 = 4H(6Hf + 0) and 1 + / — tp 7^ 6f. 
In the second case, (b), we have, 

H = ttt^' (22) 

u _ 2H(6fH - 

H -i + f-<,-er (23) 

and, 

0(4if f - j) - 4^ 

2(1 + /-^ -6/0 1 J 

This case never occurs because when H = and if 7^ we obtain -0 = Qf'H. By 
replacing it in Eq. (122j) one leads to 6/' = 1 + / — tp which contradict our primary finding. 
Finally, for the third case, (c), we have, H and H as Eqs. (|T6|) and ( TT7|) in addition to, 



(4Hf"{6H + 0) - f"0 2 



2(l + /-^-6/') ' (25) 

Therefore the conditions are (c-1) AHf"(6H + 0) 7^ /'"0 2 , (c-2) /'" = and /" 7^ 0, (c-3) 
/" = and /'" 7^ 0, (c-4) 7^ 0, in addition to the AHf'(QH 2 + 0) - /"0 2 = ^(4ff + 0) and 
1 + /-V^6f. 

With numerical calculation, as shown in Fig. 1, by appropriately choosing model pa- 
rameters we construct a cosmological model so that crossing the phantom divide occurs at 
t > and crosses —1 around this point which is supported by observations [l^]. In some 
sense this model is similar to Quintom dark energy models consisting of two quintessence 



and phantom fields 42]. 
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Fig. 1: The graph of u plotted as function of time for f e b ^ l \ f = 0.5 and 
b = 1.5. Initial values are 0(0) = 0.5 , 0(0) = -0.02, ^(0) = 0.5 , ^(0) = 0.01. 

By choosing t = to be the bouncing point, the solution for a(t) and H(t), Eq. (fl6|) . 
(see Fig. 2) provides a dynamical universe with contraction for t < 0, bouncing at t — 
and then expansion for t > 0. 




JJ 0.001 




100 163 201) 



Fig. 2: The graph of scale factor a and H, plotted as 
function of time for jo e 6 ^, = 0.5, and b = 1.5. Initial values are 
0(0) = 0.5 , 0(0) = -0.02, ^(0) = 0.5 , ip(0) = 0.01. 

By definition, for a bounce or turn-around process to occur, one must require that at the 
pivot point d = and d > around the bouncing point, while a < around the turn- around 
point. According to Eqs. and fTTUl) . for 2( - 1H fy_^ > one can get 



p e ff > 0, p e ff < (or p e ff > 0) for the bounce (or turn- around) 



(26) 



or equivalently, the EoS parameter approaches negative values (or positive values) at the 
bounce (or turn-around) point. This shows u possibly crosses over the cosmological constant 
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boundary {to = — 1), which interestingly implies the necessity to have the Quintom matter 
for the realization of the oscillating universe under modified gravity. 

In the context of a spatially flat four-dimensional background FRW metric, if we are 
to obtain a smooth transition from a contracting universe into an expanding phase, there 
must be a period when NEC is violated . In this case, we need a kind of matter which 
admits an EoS parameter which is less than -1, but only around the bounce. Neither 
regular not Phantom matter alone can achieve a transition in the EoS parameter through the 
cosmological constant boundary. Therefore, a Quintom model is the only possible solution 
to resolve this difficulty. 

A detailed examination on the necessary conditions requires for a successful bounce shows 
that during the contracting phase, the scale factor a(t) is decreasing, i.e., a < 0, and in the 
expanding phase we have a > 0. At the bouncing point, a — 0, and so around this point 
a > for a period of time. Equivalently in the bouncing cosmology the Hubble parameter 
H runs across zero from H < to H > and H = at the bouncing point. A successful 
bounce requires that the following condition should be satisfied around bouncing point, 

"=-2^ 1+W >' >0 - (27 > 
From Fig. 1 and 2, we see that at £ — )• 0, oj < — 1 and H is positive which satisfies the above 
condition. Also we see that at the bouncing point where the scale factor a(t) is not zero we 
avoid singularity faced in the standard cosmology. 



3. MODEL RECONSTRUCTION 



Reconstructing of the model for the EoS parameter and deceleration parameter in three 
forms of parametrization 43[ has been studied in this section. From effective energy density 
and effective pressure, Eqs, (fT3|) and (|T4|) one can define construction function K as: 



M 2 

ZPeff ~Peff= ; , / ~ W + 6(/ + 3Hf)) 



1 + 

With comparison to Eq. f|T3|) and using Eq. f lT6|) we have 



IK- 



Peff = K-3V, 



where, 



V 



Ml 



1 + 



f + AHf-Hi/j 



(28) 



(29) 



(30) 



One then simply finds that the energy pressure in terms of new functions K an V as, 

p eff = K - V, (31) 
now we can rewrite the EoS parameter as, 

(2K-AV\ , . 

u = -l+ ^ =- ■ 32 

\K-3V ) 

It can be seen that u > — 1 when K > 3V or K < 2V and u < — 1 when 21/ < K < 3V 
with the constrain that K ^ 3V. The transition from w > -1 to w < -1 happens when 
K = 2V. 

In order to reconstruct the cosmological parameters we rewrite the modified Friedman 
equations in the present of dust matter as, 

3M 2 p H 2 = Pm + p eff = p m + K-3V, (33) 

-2M 2 p H = Pm + Peff + Peff = Pm + 2K - 4V, (34) 

where p m is the energy density of dust matter. We thus have new K and V in the present 
of matter: 

k = ^ Pm -3M 2 (2H 2 + H), (35) 
V = ^ Pm -M 2 (3H 2 + H)- (36) 



In here, we assume that the two effective and dust matter fluids do not interact. From [44 1 
the expression for the energy density of dust matter with respect to the redshift z is given 
by, 

p m = 3M 2 H 2 n m0 (l + z)\ (37) 

where Q m o is the ratio density parameter of matter fluid and the subscript indicates the 
present value of the corresponding quantity. One then can rewrite K and V with respect to 
the redshift z as, 

K = ^M 2 H 2 (fi m0 (l + zf - 4r + (1 + z)rW) , (38) 
V = X -M 2 V E 2 (3Q m0 (l + zf - 6r + (1 + z)r^) , (39) 
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where r — ^ and = By using Eqs. ( 138]) and ( 139]) the EoS parameter can be 

rewritten as, 

(1 + z)rM - 3r 

where then r(z) can be evaluated 

r(z) = Q m0 (l + zf + (1 - n m0 )e 3 fo zl -W ld z. (41) 

Moreover, by employing 4 = — (1 + 2:)^-, the deceleration parameter g can be obtained as, 

_ * (l + z^-ir, 
yV 7 # 2 2r V ; 

Now, with the following three different forms of parametrization, using numerical calcula- 
tion, we reconstruct EoS and q(z) parameters. 

Parametrization 1: 



This parametrization has been proposed by Chevallier and Polarski 



45] and Linder 46] 



where the EoS parameter of dark energy in term of redshift z is given by, 

u(z) = ^ + ~ (43) 

By fitting this model to the observational data we find that Q m o = 0.29, w = —1.07 and 
w fl = 0.85 

Parametrization 2: 

The EoS parameter in term of redshift z has been proposed by Jassal, Bagla and Padman- 

as, 



abhan 
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u{z) = UJo+ oTW' ( } 

vhere again fitting the data [47], f2 m0 = 0.28, c^o = —1.37 and U}, = 3.39 
Parametrization 3: 



The third parametrization has proposed by Alam, Sahni and Starobinsky 49], where r(z) 
given by, 

r(z) = a rn0 (l + zf + A + A^l + z)+ A 2 (l + zf- (45) 
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This parametrization can be thought as the parametrization of r(z) instead of u(z). By 
using the results in j^], we get coefficients of the this parametrization as f2 m o = 0.30, 
A = 1, A 1 = -0.48 and A 2 = 0.25. 

The evolution of oj{z) and q(z) are plotted in Fig. 3 for the three form of parameteri- 
zations. The figure shows that the EoS parameter crosses the phantom divide line for the 
first and second parametrization and never crosses the line in the third parametrization. 
The second parametrization crosses the phantom line in two different values of z. As can be 
seen from the graph of deceleration parameter, for the second parametrization the universe 
undergoes an acceleration period from z = 0.39 until now in comparison to the first and 
third parameterizations that acceleration starts earlier. 




Fig. 3: The graphs of the EoS parameters, uj, and deceleration parameters, q, with 
respect to the redshift z. The solid, dot and dash lines represent parametrization 1, 

2 and 3 respectively. 

Also, using Eqs. (|38|) . (139]) and the three parameterizations, the evolutions of K{z) and 
V(z) are shown in Fig. 4. 
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Fig. 4: The graphs of the reconstructed K and V with respect to the redshift z. The 
solid, dot and dash lines represent parametrization 1, 2 and 3 respectively. 

Since K and V are now known functions of z, we can obtain the evolutions of <f>(z) and 
f(z) with respect of z, which are plotted in Fig. 5 for the three parameterizations. One 
may also directly obtain the relationship between the function / and the scalar field 0, 
which is plotted in Fig. 6 for the three parameterizations. 
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Fig. 5: Graphs for the reconstructed / and m respect of z. The solid, dot and 
dash lines represent parametrization 1, 2 and 3 respectively. Initial values are 
0(0) = 0.5 , 0(0) = -0.02, ?/>(0) = 0.5 , ^(0) = 0.01 and /(0) = 0.5. 

One then can reconstruct f(4>), which is plotted in Fig. 6. The exponential behavior of 
the reconstructed /(0)with respected to in particular for the second parametrization is 
compatible with the one initially assumed in the numerical calculations. 
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Fig. 6: Graphs for the reconstructed / in respect of 0. The solid, dot and 
dash lines represent parametrization 1, 2 and 3 respectively. Initial values are 
0(0) = 0.5 , 0(0) = -0.02, ^(0) = 0.5 , ^(0) = 0.01 and /(0) = 0.5. 



4. SUMMARY AND CONCLUSION 

In this paper , we consider a local scalar-tensor formulation of non-local gravity as a 
simple modified model characterized by two scalar fields and ip and function /(0) which 
can be viewed as scalar potential in the model. Analytical study of the solution shows 
that under special condition, the universe may go through a transition from quintessence to 
phantom phase which is also supported by numerical analysis. 

In analytic studying of the dynamics of the EoS parameter we achieve the constraints that 
one has to impose on the scalar fields and their first and second derivatives in order to have 
phantom crossing. In numerical approach, the EoS parameter crosses u = —1 for t > 0. 
We also find that universe undergoes a bounce, i.e contracts, reaches a minimum radius and 
then expands. 

Finally, we reconstruct our cosmological parameters, such as the EoS parameter, decel- 
eration parameter and potential function /(0). In general, different cosmological models 
can be observationally differentiate in terms of the potential function of the dynamical sys- 
tem. In here, in the construction of the deceleration parameter d(z), it is found that the 
strongest evidence of acceleration occurs about redshift z ~ 0.2-0.39 |47B. Also for the EoS 
parameter, it is found that U)(z) crossing to be around redshift z ~ 0.2 j^]. This suggest 
that all three parametrization are suitable for small z and in particular the second one is 
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remarkably close to the observational data. We then reconstructed the potential function 
of the dynamical system using three forms of parametrization by fitting the model to the 
observational. The behavior of f((p) is similar to the initially assumed exponential form. 
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